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Abstract 



In a large class of supersymmetric (SUSY) axion model the mass of axino a 
(a fermionic superpartner of the axion) is predicted as m~ < 0(1) keV. Thus, 
the axino is the lightest SUSY particle (LSP). We pointed out that such a light 
axino provides a natural solution to the gravitino problem, if the gravitino 
is the next LSP. We derive a constraint on the reheating temperature Tr of 
inflation, Tr < 10 15 GeV for the gravitino mass 1713/2 — 100 GeV, which is 
much weaker than that obtained in the minimal SUSY standard model. 
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The CP violation in QCD is one of the most serious problems in the standard model. 
In spite of continuous effort to solve the strong CP problem in the last coupled decades, 
the mechanism proposed by Peccei and Quinn |]J is still the most attractive one. The 
spontaneous breakdown of the Peccei-Quinn symmetry produces a Nambu-Goldstone boson 
(called as axion "a") || and the breaking scale F a is stringently constrained by laboratory 
experiments, astrophysics and cosmology as F a ~ 10 10 -10 12 GeV ||. 

Supersymmetric (SUSY) extension of the Peccei-Quinn mechanism necessarily predicts a 
fermionic partner of axionfjso-called axino a, whose mass is highly dependent of models 
However, in a large class of SUSY axion models || the mass of axino is predicted in the 
region m~ < 0(1) keV which is cosmological harmless [||]. In these models the axino a is 
the lightest SUSY particle (LSP) and the gravitino ^3/2 can decay into a pair of the axion a 
and the axino a. We point out, in this letter, that the light axino provides a natural solution 
to the cosmological gravitino problem [0, if the gravitino ^3/2 is the next LSP. We assume 
that the axino a is the LSP of mass m~ < 0(1) keV and the gravitino is the next LSP of 
mass m.3/2 — 10 2 GeV throughout this letter.^ 

Radiative decays of gravitinos are cosmologically dangerous, since they take place after 
the epoch of the big-ban nucleosynthesis (BBN) and destroy light nuclei synthesized by the 
BBN. To avoid this gravitino problem, the reheating temperature Tr of inflation should be 
low enough f7|. It is found in Ref. |J that the reheating temperature should be T R < 10 6 GeV 
for m 3/2 ~ 100-500 GeV and T R < 10 8 GeV for m 3/2 ~ 500 GeV-1 TeV. Such a low reheating 
temperature involves significant physical implication, i.e., it excludes some inflation models 
and/or baryogenesis scenarios. The relevant example here is the leptogenesis []TD via decays 
of heavy Majorana neutrinos to account for the baryon asymmetry of the present universe. 
For the case when heavy Majorana neutrinos are produced by thermal scatterings, which is 
the most conventional production mechanism, a successful leptogenesis requires the cosmic 
temperature of 10 10 GeV [fLljj , which leads to the gravitino problem. One of motivations 
in this letter is to solve this problem by relaxing the above constraints on the reheating 
temperatures. As shown below, our hypothesis, i.e, the axino is the LSP of m~ < 0(1) keV 
and the gravitino is the next LSP of m 3 / 2 — 10 2 GeV, allows the reheating temperature of 
10 15 GeV which is sufficiently high for the thermal leptogenesis to work.0 

In the present model, the main decay of the gravitino is ^3/2 a + a and its lifetime 
is estimated as 



1 The axion supermultiplet $ can be written by $ = a + ia + V29a + 9 2 F§, where a denotes an 
axion, a a saxion, and a an axino. 

2 This possibility was considered in the context of the galaxy formation (8|. However, the gravitino 
is assumed to have a much longer lifetime than the estimate in Eq. (||), and hence their analysis in 
not applicable for the present purpose. Furthermore, cosmological constraints on the SUSY axion 
model discussed in this letter were not investigated there. 



3 Another solution had been proposed in Ref. 12]. 
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where M* = 2.4 x 10 18 GeV is the reduced Planck mass. Since a and a have very weak 
couplings to the ordinary particles, this gravitino decay does not destroy any BBN light 
nuclei. The ratio of the gravitino energy density to the entropy density is given by 



PS/2 1n _ 9 n / m ~g V ( T R \f m 3/2 



-1 



where m~ is the gluino mass. The ratio Pz/2/ s in Eq. (|2|) should be smaller than about 10~ 4 
GeV. Otherwise, this extra energy density raises the expansion rate of the universe at the 
BBN epoch and leads to overproduction of 4 He. This gives an upper bound on the reheating 
temperature as T R < 10 15 GeV for m 3/2 = 10 2 GeV.Q 

On the other hand, the lightest SUSY particle \l next to the gravitino decays into a 
gravitino emitting photons. If this is only the decay mode, the energetic photons destroy 
the light nuclei and cause a serious problem in the BBN |7|. This is because the decay takes 
place soon after the BBN ends. However, in the present model such a particle can decay 
mainly into an axino and a photon,[|and its decay lifetime is |T5|p| 
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That is, it decays much before the BBN starts and hence there is no problem at all. 

Thus, the problem we must discuss below is whether the SUSY axion model with the 
light axino is cosmologically safe or not. First, we discuss cosmological abundance of a, 
especially, the overclosure problem of the LSP axino. 

Let us discuss possible production mechanisms of the axinos a. In the early universe the 
axinos are produced in the thermal equilibrium though the reactions like qq *-> ag, and it 
decouples from the thermal bath at the cosmic temperature |J 

T *~ lo8GeV (ira) 2 - < 4 > 

If the reheating temperature of inflation is higher than this decoupling temperature (Tr ^> 
Ta), the yield of the axino Y~ (Y~ = n~/s with the axino number density n~ and the entropy 
density s) is estimated as 



4 A similar constraint on Tr was obtained for the lighter gravitino of mass ~ 100 MeV in Ref. || 
from their scenario of the structure formation. However, our condition from Eq. (Q) leads to a 
more stringent constraint on the reheating temperature Tr < 10 12 GeV for such a light gravitino. 

5 xl is assumed to be mainly composed of the photino 7. 

6 If the -R-parity is broken, xl can decay into the ordinary light particles avoiding the problem in 
the BBN. 
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Y~a^ n ~~ 1(T 3 • (5) 



s 



On the other hand, for Tr <C T d , the yield of the axino is given by 

Y~ a ~ 10- 3 (^) . (6) 

For the case of the stable LSP axino, the present energy density of the axino may exceed the 
critical density of the present universe in some parameter regions. The density parameter 
of the axino is 

° P J So ' 

where p c is the critical density and So denotes the total entropy density of the present universe 
(Pc/ sq = 3.6 x 10~ 9 h 2 GeV with the Hubble parameter h in unit of 100 km/sec/Mpc). From 
Eq. (H) we find 

n~h* ~ 5.8 x 10 5 ( T ^ v ) . (8) 

Therefore, the non-overclosure limit fi~h 2 < 1 gives the upper bound on the axino mass 
as[§ 

m~ < 2 keV . (9) 
On the other hand, if Tr < T d , we find that the upper bound (||) is relaxed as 

m~ < 2 keV (^) . (10) 

It should be noted here that we have a large class of SUSY axion models with such a 
light axino, as mentioned in the introduction. 

The axinos are also produced by the decays of gravitino and xl- However, they are safely 
neglected because the axino mass should be small enough to satisfy the condition Eq. (||) or 

©• 

Next, we turn to discuss the cosmological problem associated with the saxion a. The 
saxions are also produced through the thermal scattering processes as well as the axinos, 
and its decoupling temperature is also given by Eq. (|]). Therefore, the yield of the saxion 
is estimated as 

/ 10 3 for T i? »T d 

n ~\l0- 3 (§) T R «T d ■ 

Then, the ratio of the saxion energy density to the entropy density is given by 

p a f 10- 3 m a for T R > T d 



s ~ \ 10~ 3 m CT (^) for T R < T d 



(12) 
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Notice that the saxion mass is comparable to the gravitino mass (m a ~ m 3/ / 2 ). For 
m a ~ 10 2 GeV the saxions dominate the energy density of the universe after the cosmic 
temperature T cools down to ~ 100 MeV. However, the saxion is not stable. The relevant 
decay channels are a — > 2g and — > 2a, and their decay rates are estimated as 

T - 2g ~ 32^if' (13) 

r CT _2a = , (14) 

where C is the constant of C ^ 1. Since the constant C depends on the model for the 
U(1)pq symmetry breaking, we take it as a free parameter. When the constant C is large 
enough as 



the saxions decays mainly into axions much before they dominate the energy of the universe 
and the produced axions are small enough. On the other hand, when C becomes smaller 
than this critical value, the a — > 2a decay channel should be suppressed enough, otherwise 
the extra energy density of the produced axions at the cosmic temperature T ~ 1 MeV 
spoils the success of the BBN. In order to avoid this difficulty, the branching ratio of the 
saxion decay into two axions should be smaller than about 0.1, i.e., C ^$ 10 -4 . In this case 
the saxions dominate the universe before they decay, the universe is reheated again by the 
saxion decay. The reheating temperature T a is estimated as 

/ m a \ 3 / 2 /10 11 GeV\ , . 

Therefore, the saxion decay completes before the BBN starts. If the Peccei-Quinn breaking 
scale is large as F a ~ 10 12 GeV, the saxion decay increases the entropy density of the universe 
by the factor A 0Q 

A~24fl^V /2 ( F - ) . (17) 
\ m a ) VlO 12 GeV/ 1 ] 

However, there is no entropy production by the saxion decay for the case of F a ~ 10 10 -10 11 
GeV. For T R <C T d , the entropy production rate A is suppressed by the factor (T R /T d ) and 
no entropy production takes place when T R /T d < 0.04. 

Furthermore, it should be noted that the saxion may be produced effectively in the 
form of the coherent oscillation after the inflation. We assume here that the supergravity 
effects induce positive mass squared for Peccei-Quinn scalar fieldsQof order of Hi during the 



If there is an entropy production after the QCD phase transition, the upper bound on F a is 
raised up above F a ~ 10 12 GeV jOJ. 

8 Peccei-Quinn fields are scalar fields responsible for the Peccei-Quinn symmetry breaking. 
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inflation (Hj denotes the Hubble parameter for the inflation). Thus, it is quite natural that 
the Peccei-Quinn symmetry is restored during the inflation if Hj > T a |]The Peccei-Quinn 
symmetry breaking occurs when the Hubble parameter becomes comparable to the scale F a , 
and the coherent oscillation starts. Then, the ratio of the energy density of the coherent 
oscillation p osc to entropy density for T <C Tr is estimated as 

Oosr- 1 Tf>F^ 



Tr W F a 



s 8 Mq ' 

Comparing it with the energy density of the saxion produced by the thermal scatterings in 
Eq. (0), we can safely neglect the energy density of the oscillation if Tr < 10 15 GeV for 
F a ~ 10 11 GeV. 

Finally, we should comment on the axion domain walls. We have assumed that the 
U(1)pq symmetry is restored during the inflation, and hence the axion domain walls might 
be formed after the inflation ends. However, this can be easily evaded by adopting a hadronic 
axion model [[T7[] with the domain wall number N^w — 10- 

In this letter we have pointed out that the cosmological gravitino problem can be solved 
by the SUSY axion model which is the most natural solution to the strong CP problem. In 
the present model, the axino is the LSP and the gravitino is the next LSP. The gravitino 
decays into a pair of the axino and the axion eluding the photo-dissociation constraint on the 
reheating temperature Tr < 10 6 -10 8 GeV which was obtained in the minimal SUSY standard 
model with an unstable gravitino of 7773/2 — 100 GeV-1 TeV We find that the reheating 
temperature can be high as Tr ~ 10 15 GeV for 7713/2 — 10 2 GeV. Therefore, the present model 



makes thermal leptogenesis scenarios fl 1] to work well without the cosmological gravitino 
problem.^ 
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Since the U(1)pq symmetry is restored during the inflation, there is no massless mode and hence 



no isocurvature fluctuation is generated [16|. 



10 The primordial lepton (baryon) asymmetry is not diluted away by the saxion decay, since we 
have only a small entropy production as shown in Eq. (|l7|). 
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